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ABSTRACT 

We argue that X-ray flaring variability observed in the transient X-ray pulsar A0535-I-26 is due to low-mode magnetospheric in- 
stability. This instability develops at the onset of accretion, in the thin boundary layer between the accretion disk and neutron star 
magnetosphere. As a result, the matter collected in the boundary layer can rapidly fall onto the NS surface close to the magnetic 
poles, but not exactly along the field lines by which the stationary accretion proceeds. This explains the shift in cyclotron line energy 
measured using RXTE data in a pre-outburst spike, with respect to the line energy observed during the main outburst. Furthermore, 
the instability can account for the difference in pulse profiles, and their energy evolution that is dilferent in the pre-outburst flare and 
main outburst. 

Key words, accretion, accretion disks - stars: neutron - X-rays: binaries 



1. Introduction 

The transient X-ray binary pulsar A0535-)-26 was discovered 
by the Ariel V satellite during a giant outburst jCoe et al. 19751 
[Rosenberg et al. 1975| l. The neutron star (NS) rotates with a 
period P ^ 103 s and is in an eccentric (e - 0.47) or- 
bit about a massive 09.7IIIe star with an orbital period of 
Porb ~ 110 d (ICoe et al. 2006t . The X-ray activity of this 
transient source is complicated: several giant outbursts were 
observed in 1975, 1980, 1983, 1989, 1994 and 2005, with 
weaker outbursts sometimes observed at successive perias- 
tron passages. This type of X-ray activity is common for 
Be/X-ray binary systems, which represent the largest subclass 
of massive X-ray binaries. The giant (type II, according to 
IStella et al. 1986t outbursts can take place at any orbital phase 
and most probably are triggered by the enhanced activity of the 
optical star dCoe et al. 20061 1. More regular, weaker type I out- 
bursts are associated with enhanced interaction of the neutron 
star with a circumstellar disk surrounding the Be-companion 
( [Okazaki & Negueruela, 200 l| l. A0535-I-26 did not show out- 
bursting activity between 1994 and 2005. The last giant out- 
burst of A0535+26 occurred after a prolonged quiescence pe- 
riod in May-June 2005 (ITueller et al. 20051 1. In September 2005, 
the subsequent normal type I outburst was extensively studied 
by RXTE and INTEGRAL (ICaballero et al. 20071 2008). These 
observations revealed somewhat unexpected features, which we 
summarize below. 

Feature 1. A short (fraction of a day) pre-outburst X-ray 
spike, with a luminosity comparable to the outburst maximum, 
is observed during the rise of the outburst. 

Feature 2. The NS pulse period at the onset of the outburst, 
remained constant within the margins of error, but rapidly de- 
creased after the periastron passage. 
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Feature 3. RXTE data shows that the energy of the cy- 
clotron resonance scattering feature (CRSF) during the spike at 
90% confidence is 52.0+| ^ keV Caballero et al. 2008l l. which 
is notably higher than the value 46. 1 + 0.5 keV measured dur- 
ing the main outburst at the same level of X-ray luminosity. 
The latter value is consistent with INTEGRAL measurements 
(ICaballero et al. 2007b . 

Feature 4. The CRSF energy in the main outburst (within the 
margins of error) is independent of the X-ray luminosity. 

Feature 5. The X-ray pulse profiles during the spike are dif- 
ferent from the pulse profiles observed during the main outburst. 
Pulse profiles during the outburst show an abrupt change at en- 
ergies above the cyclotron resonance, while those in the pre- 
outburst spike do not. 

The additional inspection of the Swift BAT (15-50 keV) light 
curve of the source (Fig. 1) reveals that the pre-outburst spike 
observed by RXTE is only one of a collection of flares, with the 
characteristic time up to a few times lO'* s, appearing at the ris- 
ing part of the outburst. The flux evolution then proceeds more 
smoothly during the spin-up of the NS {Feature 6), which be- 
gins close to periastron passage. A similar flaring activity was 
observed during the onset of the preceding giant outburst, at the 
same flux level, in April 2005, and at the next periastron passage 
in December 2005. 

In this Letter we argue that the flaring activity and pre- 
outburst spikes are most likely due to magnetospheric instabil- 
ity. This would take place at the onset of a smooth increase in 
the mass-accretion rate through the accretion disk surrounding 
the NS, close to periastron passage. This instability would cause 
plasma that had accumulated in the thin boundary layer, between 
the accretion disk and NS magnetosphere, to fall onto the NS 
surface as large blobs. These blobs would be channeled to re- 
gions close to the NS magnetic poles by different magnetic field 
lines than those guiding the accretion flow smoothly increasing 
from the disk. 
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Fig. 1. The Swift/BAT light curve of A0535+26 during the normal outburst in August/September 2005. The smooth solid line shows 
the pulse period development as measured by RXTE (from Caballero et al. 2007b). The dashed line shows the time of periastron 
passage. 



2. The model 

Since the 1970s, it has been known that matter can penetrate 
the NS magnetosphere, and be accreted onto magnetized NS, 
via various instabilities. The gravitational interchange instabil- 
ity (Rayleigh-Taylor, or Kruskal-Schwarzschild [KS] for plasma 
in magnetic field) of cold plasma at the magnetospheric bound- 
ary (see lElsner & Lamb 19761 lArons & Lea 1976al l. in partic- 
ular, mediates quasi-spherical accretion. In this case, close 
to the magnetic equator where KS instability is more likely, 
plasma enters the magnetosphere as bubbles or filaments that 
fall deeper towards the NS. All or part of these filaments, bro- 
ken up into smaller pieces by Kelvin-Helmholtz instabilities 
( lArons & Lea 1976al l, can then be entrained by the magnetic 
field and channeled toward the magnetic polar region on the NS 
surface. 

During disk accretion onto a rotating magnetized NS the 
Kelvin-Helmholtz instability between the disk material and the 
magnetosphere is also important dAnzer & Borner 19801 1983). 
The transition between the accretion disk and the rotating NS 
magnetosphere (for the aligned magnetic dipoles) was studied 
by Ghosh & Lamb (1978, 1979), Scharlemann (1978), Anzer & 
Borner (1983) and later by Lovelace et al. (1995) using more 
realistic assumptions. In the latter model, the twisting of the 
stellar magnetic field by the differentially rotating disk leads 
to the appearance of open magnetic field lines extending out- 
ward from both the star and the disk. The accretion process 
may then be maintained by an MHD outflow from the disk that 
transports away angular momentum. Numerical MHD simula- 
tions of the disk interaction with an inclined magnetic dipole 
( IRomanova et al. 2"0 03) confirmed the basic features of the disk- 
magnetospheric interaction. Accretion onto NSs can be unsta- 
ble, leading to outbursts and flares, because of instabilities in the 
magnetosphere (IBaan 19791 Spruit & Taam 199311. The KS in- 



stability affects accretion onto magnetized NSs at both large and 
small accretion rates. At small accretion rates, matter accumu- 
lates at the boundary layer and flaring events are expected. At 
large accretion rates, matter is accreted more steadily. 

In the case of A0535+26 we deal with non-stationary accre- 
tion triggered by enhanced NS interaction with the circumstellar 
disk and the wind of the Be-companion at binary periastron pas- 
sage. BATSE observations of this source JBildsten et al. 1997] l 
indicated the strong disk accretion spin-up during type II (giant) 
outbursts, while no significant spin-up was observed during a 
series of type I outbursts, at successive periastron passages both 
before and after the giant outburst. Moreover, the average spin- 



down of the NS rotation is clearly seen in the BATSE data be- 
fore and after the giant outburst, suggesting that a propeller-like 
mechanism operates between accretion episodes close to perias- 
tron. We expect similar long-term spin-up/spin-down behavior 
of magnetized NSs in eccentric Be X-ray binaries on very gen- 
eral grouds JStella et al. 19861 1. The magnetic propeller regime 
was investigated by Lovelace et al. (1999). Numerical MHD 



simulations (Ustyugova et al. 2006 1 revealed the non-stationary 



character of mass outflow during the propeller stage. Clearly, 
the transition from propeller to accretion stage is the most com- 
plicated for the analysis. Here we suggest that the non-stationary 
features observed in the outburst of A0535+26 relate to the mag- 
netospheric instability which develops during the transition to 
the accretion stage. 

The pre-outburst X-ray flaring activity is not an exceptional 
feature of the 2005 outbursts in A0535+26. Hints of other flares 
in A0535+26 are apparent for several other type 1 outbursts in 
BATSE light-curve data for its giant outburst in 1994 (see Fig. 
1 of Finger et al. (1996)). Similar pre-outburst spikes are ob- 
served in other Be/X-ray binaries, e.g. in BATSE observations 
of GS 1843-02 ( [Finger et al.T999| and INTEGRAL observa- 
tions of EXO 2030-H375 dCamero Arranz et al. 2005T l. A plau- 
sible picture may be delineated as follows. 

1) We assume that an accretion disk surrounds the NS in 
A0535+26. It is definitely the case during the giant outbursts, 
as the spin-up measurements and QPO observations suggest 



(Sembayetal. 1990 Finger etal. 1996 1. The disk can survive 



for several orbits after the giant outburst. The observed strong 
spin-up during the main part of the August/September 2005 out- 
burst (Caballero et al. 2008 and Fig. 1) confirms the presence of 
the disk. 

2) When the accretion rate onto the NS increases as it ap- 
proaches orbital periastron, the magnetospheric radius R,,,, de- 
termined by magnetic pressure balancing plasma pressure, de- 
creases according to the proportionality/?,,, oc M^'^l'' . Accretion 
occurs when the stopping radius of matter Rg, generally of 
the order of becomes smaller than the corotation radius 
Rc = {GMIap-yi^, where M is the NS mass and w = InlP is 
its spin frequency. In all models of the disk-magnetospheric in- 
teraction a thin boundary layer between the disk and the rotat- 
ing NS magnetosphere does exist in which the matter becomes 
gradually captured by the magnetic field lines. The width of this 
region is model dependent, but is of the order of the disk thick- 
ness (e.g., < IQcs/cjk in lAnzer & Borner 19831 or ~ Scs/cok 
in ILovelace et al. 19951 here Cj is the sound speed and otK is 
the Keplerian frequency). So for an estimate we can assume 
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Al ~ h ^ O.lRa- The amount of mass accumulated in this 
layer can be roughly evaluated assuming the standard o'-theory 
([Shakura & Sunyaev 1973]l: 
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AM = p2nR„2hAl - (4 x lO'^g)Rl[^a-^''Ml[l{l - ^ 

where /?„ 9 is the stopping radius in units of 10"^ cm, M-\q is the 
mass accretion rate in units of 10"'° Mg/yr, ^ - J KM y/GMRa) 
is the angular mometum flux through the disk normalized to the 
flux from the stopping radius. Here we have assumed the gas 
pressure to dominate over the radiation pressure in the boundary 
layer and the main opacity to be due to Thomson scattering. Of 
course, the disk is already not Keplerian at this radius [doj/dr — » 
at Ra), but for the case of interest Ra ~ Rc we neglect the 
deviation from the Keplerian motion. At the onset of accretion 
^ - 0, the viscosity parameter can be very small, a ~ 0.01 or 
smaller, so the amount of matter in the boundary layer may be 
as high as ~ 10^' g. 

3) As inferred from the analysis by Baan (1979), at small 
accretion rates (~ lO""^M0/yr = lO'* g/s) the magnetopause in 
A0535+26 with P ~ 104 s may be on the verge of a KS instabil- 
ity. The model by Baan (1977, 1979) has sucessfully described 
the statistical properties and the burst waiting time - luminos- 
ity correlation of type II outbursts in the rapid burster MXB 
1730-335. In the linear regime, the KS instability grows with 
time exponentially in proportion to exp(rf) with the increment 

= 't^etf tanh(fe) dArons & Lea 1976bl|Baan 1979D , where k is 
the angular wave number, z is the characteristic height above the 
interface and ^eff is the effective gravity acceleration, geif takes 
into account the gravitational attraction (oc l/R^,), minus sta- 
bilizing contributions from the centrifugal force (oc ci/R„), the 
bending of field lines by the disk, the curvature of the magnetic 
field near the magnetopause (oc 1 /(pmRl,)), possible velocity gra- 
dients in the plasma flow, etc. The condition for the instability to 
occur is the change of the sign of ^eif from negative (stability) to 
positive (instability), which defines the zero-surface ^eff (^;) = 
at some distance from the NS. We note that under the assumption 
of the corotation of matter and the magnetosphere, the instabil- 
ity can occur only when the zero-surface is inside the corotation 
radius, R~ < R^ because of the centrifugal term in ge/y- This 
implies that instability can develop only when accretion is per- 
mitted centrifugally. This is consistent with the observation in 
A0535+26 of flaring activity, and accretion. 

As the magnetopause moves inwards through the zero- 
surface with increasing M, the simplified analytical treatment 
by Baan (1979) shows that the lowest harmonics are the first to 
become unstable. This finding however remains controversial. 
Early numerical simulations of the disk-magnetosphere interac- 
tion (e.g. Rastatter & Schindler 1999) have shown that the high- 
m modes grow faster. The contribution of the component of 
the field in the disk can suppress high-m modes. Simulations by 
Wang & Nepveu (1983) have shown that the high-m modes grow 
faster but may later merge, forming larger filaments. The first 
global 3D MHD simulations of the disk-magnetosphere inter- 
action through the KS instability (Romanova et al. 2007) have 
shown that the low-m modes dominate for a wide range of pa- 
rameters. It is not however clear which factors are responsi- 
ble for suppressing the high-m modes. These 3D simulations 
have shown that at small accretion rates matter accretes through 
funnel streams, while at large accretion rates, matter accretes 
through gravitational KS instabilities. The 3D simulations have 
been carried out for quite small magnetospheres (~ IORns) 
which are more applicable for the millisecond pulsars. This is 
probably why the accumulation of matter at small M and burst- 
ing events have not yet been observed in these simulations. 



Therefore, we can assume that for Be/X-ray transients with large 
magnetospheres, such as A0535+26, low-mode magnetospheric 
instability, initiated by an increase in M, can develop. 

In such binary systems the infall of plasma blobs from the 
magnetopause can also be triggered by variation of the rate 
of mass captured from the variable stellar wind. So short X- 
ray flares or spikes can in principle be observed at any orbital 
phase. An "unexpected" X-ray flare in A0535+25 observed by 
INTEGRAL in October 2003 during the prolonged quiescence 
period ( |Hill et al., 2007| i could be an example of such an event. 

4) The low-mode magnetospheric instability can break-up 
the entire boundary layer in the disk. All material stored in the 
boundary layer may then rapidly enter the magnetosphere. Our 
estimate above indicates that up to ~ 10^' g can fall onto the NS 
surface on the free-fall time scale (of the order of 20 s) producing 
short X-ray spikes with a maximum X-ray luminosity of up to 
~ 5 X lO-'^ erg/s (feature 1 from the Introduction). The amount of 
matter falling onto the NS surface during the spikes in the con- 
sidered outburst of A0535+26 can be estimated from the light 
curve (Fig. 1) and is found to be AM ^ j Mdt < 10^' g in each 
spike. The characteristic timescale of the spikes is ~ 10'' s and 
agrees with the time required to replenish the boundary layer 
~ AM/M. The angular momentum supply from the disk dur- 
ing flares changes the NS spin period by the fractional amount 
|AP/P| = AM \1GMR„PI{2kI) < 7 x lO"*" where we assume that 
the NS moment of inertia is / = lO''^ g cm- and /?„ = /?e = 10^ 
cm. This is within the uncertainty of the period measurement 
(AP/P ~ 5x lO ^ ICaballero et al. 2008b . which explains the lack 
of noticeable spin-up of the NS at the initial phase of the outburst 
(feature 2 from the Introduction); see also the BATSE observa- 
tions of A0 535-H26 presented i n jBildstenet al. 1997| i, and of GS 
1843-02 in|Fingeretal. 1999| 



5) The plasma entering the magnetosphere via the low-mode 
KS instability may become frozen into the magnetic field more 
close to the NS than the main accretion flux and hence fall 
along different magnetic field lines than those by which the 
quasi-stationary accretion is channeled. This explains the ob- 
served difference in the cyclotron line energy during the ini- 
tial X-ray spike and the remainder of the outburst (feature 3 
from the Introduction). Similarly to Her X-1 (Staubert et al., 
2007), the absence of a radiation-dominated accretion column 
in A0535-I-26 is suggested by the independence of the CRSF 
energy late in the outburst from the observed X-ray luminos- 
ity (feature 4 from the Introduction). The height of the emission 
region above the NS surface is then about several hundred me- 
ters. To change the CRSF energy by the fractional amount of 
AEc/Ec ~ 10%, the emission region during the spike is required 
to move towrds the NS surface by the amount AR/R^s ~ 3% 
(assuming the dipole magnetic field), i.e. by about 300 meters. 
Therefore, the emission from the spike is likely to originate very 
close to the NS surface. This explains the different pulse profiles 
during the spike, especially in hard X-rays (feature 5 from the 
Introduction). 

The smooth change in pulse profile shape at energies above 
the cyclotron resonance (as in Her X-1, Klochkov et al. 2008) 
may correspond to a pencil-beam diagram of emission during 
the flare. The photon cross-section may then change smoothly, 
at the cyclotron resonance energy, in the strong magnetic field 
along the direction of the field. In the main outburst, however, 
the accretion column is higher (see above), the density increases 
such that it can become optically thick in the direction normal 
to the field, and if there is a temperature gradient to the cen- 
ter, a fan-like beam can be additionally formed by extraordinary 
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photonfl The cross-sections of extraordinary photons far from 
the resonance (E <s: Ec) are proportional to ctj^ ^ crj(E /E^)^, 
which is a far smaller cross-section than that of ordinary pho- 
tons cr|| ^ cr7-(sin^0 + cos~ 9{E / Ec}^). The e-photons are then 
able to escape effectively from larger optical depth. (Here ctj- is 
the Thomson scattering cross-section and 9 is the angle between 
the incident photon and the magnetic field; see Harding & Lai, 
20043). Above cyclotron resonance, cr^ - cr|| ^ crx and photons 
escape from small optical depths. This can explain the disap- 
pearance of the large hump at pulse phase ~ 0.7 after crossing 
the CRSF energy (Fig. 3 in Caballero et al. 2008). 

6) As the accretion rate increases, the NS starts spinning- 
up and the angular momentum flux through the inner disk 
rapidly approaches its maximum possible value ~ M y/GMRa 
(see Caballero et al. 2008 and Fig. 1). The parameter ^ in the 
expression for AM above increases from zero to almost one, 
which strongly reduces the mass in the disk boundary layer and 
hence the mass available for unstable accretion. The amplitude 
of spikes rapidly decreases and we do not see strong flaring ac- 
tivity during the spin-up stage (feature 6 from the Introduction); 
see also Figs. 2 and 10 in Finger et al. 1999 for the case of GS 
1843-02. This consideration suggests that flares can reappear af- 
ter the spin-up stopping, for which there is indeed a hint in the 
light curve of the A0535+26 outburst during September 13-15, 
2005. At this time the NS had stopped spinning-up (see Fig. 1). 



3. Discussion 

The magnetospheric instability model proposed for the pre- 
outburst flares in A0535+26 is generic and can be applied to 
other transients. The prerequisite, however, is that the source 
must be on the verge of low-mode instability, which depends on 
the NS magnetic field, the spin period, the accretion rate and 
possibly other parameters (e.g. the misalignment of the mag- 
netic dipole and/or NS spin axis relative to the orbital angular 
momentum). It is possible that the transition from propeller to 
accretion stage can occur without a strong KS instability being 
observed. The model predicts changes in pulse profiles and in 
the cyclotron line energy during the short flares, as observed in 
A0535-I-26 dCaballero et al. 2008l l, which can be compared with 
observations of other sources. 

Are there other possible explanations for the observed pre- 
outburst X-ray spikes? 

The SPH modeling of an accretion disk surrounding the 
NS (Hayasaki & Okazaki, 2006) reproduces normal outbursts 
at successive periastron passages. In some instances the mod- 
eling shows a single peak preceding outburst maximum. The 
accretion disk is formed from the Roche lobe overflow of the 
coplanar circumstellar disk surrounding the Be companion star 
The important feature of these SPH simulations is the transient, 
one-armed spiral structure of the accretion disk, induced by a 
phase-dependent mass accretion. Such accretion disks can be 
responsible for mass transfer enhancement close to periastron. 
The model by Hayasaki & Okazaki, however, ignores the disk- 
magnetosphere interaction critical to transient accretion onto 
magnetized NS, and can not reproduce the observed flaring ac- 
tivity. 



' Ordinary polarized photons have electric field vector E lying in the 
plane formed by the magnetic field and the wave vector of the photon; 
the E-vector of extraordinary photons is perpendicular to this plane. 

^ A similar dependence holds for free-free absorption cross-sections 
(Kaminker et al. 1983). 



Spruit & Taam (1993) also found a viscous, disk- 
magnetosphere instability that was associated with the magne- 
tospheric boundary around the corotation radius. This instability 
results in a cyclic enhancement of the mass accretion rate on 
the viscous time scale at the magnetospheric boundary. While 
this timescale can be as short as 10"* s for the corotation radius 
10^ cm, the chaotic behaviour of flares observed in A0535+26 
and, most importantly, other features (different pulse profiles and 
CRSF energy, and absence of the flaring during the NS spin-up 
phase) imply that this model can not provide the correct expla- 
nation. 

Our model for short flares observed at the onset of outbursts 
in the transient X-ray pulsar A0535+26 is based on magneto- 
spheric instability. It succesfully explains all features of the out- 
burst observed during August-September 2005 by RXTE and 
INTEGRAL, and makes clear predictions for future observa- 
tions. The present-day accuracy of cyclotron line measurements 
in accreting neutron stars, precise timing analysis and evolution 
of X-ray pulse shapes with luminosity at different energies, have 
became the working tools to probe the non-stationary accretion 
onto magnetized neutron stars. 
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